A minimum single-band model for low-energy excitations in superconducting K^Fe 2 Se 2 
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We propose a minimum single-band model for the newly discovered iron-based superconducting KvFeoSei. 
Our model is found to be numerically consistent with the five-orbital model at low energies. Based on our 
model and the random phase approximation, we study the spin fluctuation and the pairing symmetry of super- 
conducting gap function. The (n/2, n/2) spin excitation and the rf v 2_ v 2 pairing symmetry are revealed. All of the 
results can well be understood in terms of the interplay between the Fermi surface topology and the local spin 
interaction, providing a sound picture to explain why the superconducting transition temperature is as high as to 
be comparable to those in pnictides and some cuprates. A common origin of superconductivity is elucidated for 
this compound and other high-T c materials. 

PACS numbers: 74.70.-b, 74.20.Rp, 74.25.Ha 
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The iron-pnictide materials have attracted much attention 
since their discovery [Q]] due to the high superconducting 
(SC) transition temperature and qualitatively similar phase di- 
agram with the cuprates. The ground state of the parent com- 
pound exhibits a collinear antiferromagnetic (AF) order and 
the superconductivity emerges when the magnetic order is 
suppressed with the doping. Thus the interplay between the 
magnetism and superconductivity is likely rather essential for 
exploring the mechanism of the superconductivity in this class 
of materials [see e.g., Refs. J^Ht]]. The magnetism in the par- 
ent compound can be understood from an inter-band Fermi 
surface (FS) nesting picture, namely, the nesting between the 
hole and electron bands could give rise to the collinear AF 
magnetic order |4j,l5]. On the other hand, the collinear AF or- 
der may also be obtained based on a J\ - Ji Heisenberg model 
with Ji > J\/2 [6-8]. For the doped systems, the supercon- 
ductivity with s ± pairing symmetry was proposed [9], which 
is based on a scenario that the superconductivity is meditated 
by the magnetic fluctuation and the inter-band FS nesting still 
plays a key role JSJLKj. 

Recently, a new type of iron-based superconducting com- 
pounds K A Fe2Se2 was reported, with the maximum SC transi- 
tion temperature above 30K Blll4l3ll . The crystal structure is 
similar to that of the previous 122-type iron pnictides, with in- 
tercalating potassium atoms between quasi-two-dimensional 
conducting Fe-Se layers. It is a heavily over-doped compound 
with the electron doping about 0.5 with respect to the par- 
ent compound. Both the angle-resolved photoemission spec- 
troscopy (ARPES) experiments Ill4l4l8l] an d the local density 
approximation (LDA) calculations il9H2lll indicated that the 
hole FS pockets disappear and only electron FS pockets exist. 
This seems consistent with previous ARPES experimental re- 
sults for iron pnictides that the hole-like FS pockets disappear 
at the doping 6 ~ 0.15 12211 . For the doping 6 ~ 0.5, the hole- 
like band should be well below the Fermi energy and make 
no contribution to the superconductivity. Obviously, the inter- 
band FS nesting picture, which seems to work well for iron 
pnictides, is unable to straightforwardly account for the super- 
conductivity of this compound, and thus the previous physical 
picture for iron-based SC materials meets a serious challenge. 



A second challenge is to answer why the superconductivity 
and magnetism could survive at so heavily overdoped com- 
pound. For this, it is natural and important to ask whether the 
interplay between magnetism and superconductivity is still es- 
sential for the superconductivity of this compound. Moreover, 
The pairing symmetry is also an important and critical issue to 
be addressed, which may be different from that of iron pnic- 
tides due to the absence of the hole FS pockets. Actually, 
the ARPES experiments lfl4l - [l7tl indicated that the gap func- 
tion is fully gapped for the pockets around M-point while it 
is small near the T point ffl fl5Tl . which is sharp contrast to 
the s x i v 2 -wave symmetry. Up to now, the pairing symmetry 
is still unclear: both s-wave and of- wave symmetry have been 
proposed iHHH. 



From the band calculations for this compound H19h21H . it 
seems that all of the five orbitals hybridize strongly, and thus 
it was proposed that all five d-orbitals should be considered to 
construct a model l23l l24ll . Note that the five-orbital model 
is one of the most frequently-used model to describe the iron- 
pnictide compound B29I1 . However, this model includes some 
non-essential parts that may be redundant for superconduc- 
tivity (especially those bands far away from the Fermi sur- 
face), and even make it quite difficult to accurately analyze 
and understand certain essential physics behind because too 
many unknown parameters are involved. For example, if the 
electron correlation needs to be included, both inter-orbital 
and intra-orbital interactions should be taken into account; 
while the corresponding interaction strengths are difficult to 
be determined (or adjusted) in order to figure out the rele- 
vant physics. Therefore, it is highly-demanded and desirable 
to establish a rather simplified minimum model that is able 
to capture the essential features of this compound. Motivated 
by this, a minimum two-band model was proposed for iron 
pnictides ll4l, [3"oll . and many physical properties have been un- 
derstood based on this type of two-band model, especially 
for those related closely to the low energy excitations. When 
studying the K v Fe2Se2 materials, a similar two-orbital model 
was put forward, while the obtained FS size is much larger 
than that obtained from the band calculation 13111 - Intrigu- 
ingly, as discussed below, we can refine the minimum model 
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of this compound to have one band to capture the essential 
physics with better results. As seen from the band calcula- 
tions il9l - l2lll there exist two electron-like FS sheets around 
M = (n, n) point. Around the F = (0, 0), a small FS pocket 
may exist, but it is k z dependent and electron-like. The FSs 
from ARPES experiments are consistent with the LDA calcu- 
lations while the electron pockets around the F point have a 
very low spectral weight. So it is reasonable to believe that 
only the FS pockets around M point are essential in the mini- 
mum model, while all other bands that do not cross the Fermi 
energy may be neglected. Taking into account that one unit 
cell consists of two iron irons plus the band folding effect, 
only one kind of FS pocket around the X or its symmetric 
points is actually relevant in the unfolded BZ. 

In this Letter, motivated by the above considerations, we 
propose a single-band tight-banding model as a minimum one 
for the superconducting K A Fe2Se2 [32]. We show that our 
model could have a similar band dispersion crossing the FS 
with that of the five-orbital model fitted from the LDA band 
calculations for the K c Fe2Se2 materials [23]. Thus the present 
model serves as an effective one for describing the low energy 
physics of the K x Fe2Se2. The spin susceptibility is calculated 
and analyzed based on the FS topology. It is found that the FS 
topology may connect to the J1—J2- J3 couplings to provide 
a coherent picture for the spin-density-wave (SDW) instabil- 
ity of parent compound and the magnetic fluctuation in the 
SC state. Our numerical results present a natural explanation 
for the high SC transition temperature of this material. The 
SC gap is also calculated self-consistently and the robust d- 
wave SC pairing is revealed. The af-wave symmetry can be 
explained well based on the spin fluctuation picture and the 
fermiology theory. 

We start from a t — J type model including the tight-banding 
term and the local spin interaction, which reads 
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where is taken phenomenologically as the single band 
tight-banding approximation: = -2t{coskx + cos ky) - 
At' cos k x cos k y - with t = 0.08 eV and f = -0.2 eV. The 
chemical potential fi is controlled by the electron doping 5. 
Hj is the local spin interaction. From the first principle calcu- 
lation in Ref II20I1 . we assume the spin interaction to include 
the nearest-neighbor, the next-nearest-neighbor, and the next- 
next-nearest-neighbor ones: 
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Here (ij), (ij)', (ij)" represent the summation over the near- 
est, next-nearest, and next-next-nearest neighbors, respec- 
tively. We set 7, : J 2 : h = -1 : 1.676 : 0.997 with J x < 0, 
from the first principle calculation [20]. 

The bare spin susceptibility from the tight-banding part can 
be calculated as, 
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FIG. 1 : (Color online) The band structure along different momentum 
cut. The solid line is the calculated band structure from our single- 
band model. The dashed lines are the replot of the band structure for 
the five-orbital model with the parameters in Ref. |23]. 



where f(x) is the Fermi distribution function. 

The correction of the spin fluctuation from the spin-spin 
interaction Hj is included in the random phase approximation 
(RPA), 



*(q> ^) = 
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where J q = -7i(cos k x +cos k y )-2J2Cos k x cos k y —Ji(cos2k x + 
cos 2ky), is the Fourier factor of the spin coupling term Hj. 

Considering that the SC pairing is meditated by the spin 
fluctuation, we may write the linearized eliashberg's equation, 



iA(k) = - JV(k-k') 



tanh(/te k </2) 
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with yS = l/T, We consider the spin-fluctuation as the effec- 
tive pairing potential V(q) = J^xiq, 0). Since we here address 
the pairing symmetry, we focus only on the zero-energy spin 
susceptibility, which should produce qualitatively correct re- 
sults for the pairing symmetry as usual I33LI3411 . 

We first depict in Fig. 1 the band dispersion calculated from 
the five-band model for KFe2Se2 proposed in Ref [23]. As is 
seen, there are three bands below the Fermi energy and one 
band above it, with the energies less than -0.5 eV and larger 
than 0.25 eV. The SC gap is ~ 0.01 eV from the ARPES ex- 
periments As a result, these bands (not crossing the 
FS) should contribute little to the superconductivity. It is rea- 
sonable to concentrate on the band that crosses the FS, which 
plays an essential role in the superconductivity. Considering 
that the whole electron filling is 6 + 6, the three filled bands 
contribute six electrons per site and thus the electron filling for 
the band that crosses the Fermi energy is 6. These features are 
crucially captured by the the present single-band model. The 
corresponding bare band dispersion is plotted in Fig.l, which 
almost coincide with the crossing Fermi energy band. This 
similarity suggests that our model could produce the quali- 
tatively same results as those obtained from the five-orbital 
model at low energies. 

We now study the spin fluctuation based on the above 
single-band model at the doping 6 = 0.4. As mentioned 



3 




q kjn 



FIG. 2: (Color online) (a) The bare and renormalized zero energy 
spin susceptibility with J\ = -0.1 eV along different momentum cut. 

(b) The Fourier factor J q as a function of the momentum q. Panels 

(c) and (d) are the intensity plot of the bare and renormalized spin 
susceptibility, respectively. 



above, the spin susceptibility are mainly determined from the 
two factors: one is the bare spin susceptibility xo evaluated 
from the tight-banding term, the other the RPA factor 1 - JqXo 
closely related to the local spin interaction. The behavior of 
the bare spin susceptibility is usually determined by the FS 
topology, with the maximum spin excitation occurring at near 
the wave vector connecting different FS sheets. For the renor- 
malized one, the pattern of J q , representing spin interaction 
strength in the Fourier space, plays an important role. Obvi- 
ously, a larger J q will lead to a larger spin susceptibility (The 
value of J\ in J q is usually determined from the SDW instabil- 
ity point). Since our main conclusions do not depend on 
we take Ji = -0.1 eV for illustration. The bare and renor- 
malized spin susceptibility along the two-dimensional cut are 
plotted in Fig. 2(a). The Fourier factor of the spin interaction 
J q is displayed in Fig. 2(b). Figs. 2(c) and 2(d) are the inten- 
sity plots of the bare and renormalized spin susceptibility in 
the whole BZ, respectively. As is seen, the largest spin sus- 
ceptibility occurs at the wave vector (7r/2,7r/2). This result, 
elaborated below to be consistent with the FS topology, may 
be checked by the future neutron scattering experiments. 

An important feature can be drawn from the results of the 
spin fluctuation. Notably, there exists an implicit relationship 
between the bare spin susceptibility and the term J q when 
comparing Figs. 2(a) with 2(b), namely, J q and Xo nave a 
qualitatively same momentum dependence. As a result, the 
RPA correction increases the spin fluctuation, strengthens the 
peaks, and suppresses the dips of the bare one. Since the su- 
perconductivity is expected to be meditated by the spin fluc- 
tuation, a stronger fluctuation means generally the stronger 
superconductivity. The matching of the vertex term J q and 
the bare band spin excitation provides a natural explanation 
for why the SC transition temperature is high for this com- 
pound. On the other hand, due to this implicit relationship, 
the understandings of the appearance of SDW (AF) instabil- 
ity based respectively on the fermiology and the local spin 
interaction picture are also consistent. Actually, for the pre- 
vious high-T f SC materials, a similar relationship exists, e.g., 
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FIG. 3: (a) The FS in the whole BZ. (b) The replot of FSs in the 
upper right 1/4-BZ with doping 6 = 0.4 and 0.5, respectively. The 
dashed lines are the (n, 0)-shift of the FSs from the upper left 1/4-BZ 



the bare spin excitation is peaked at {n, tt) for cuprates and 
consistent with the factor J q = -J(cosq x + cosg,) when 
the nearest-neighbor spin interaction is considered. For iron- 
pnictide, J q in the J\ - J2 model is usually written as J q = 
-Ji(cosq x + cos^y) - 2J2COsq x cosq y , which has the maxi- 
mum at (tt, 0) or (0,7r) if J2 > J\/2, notingthat J2 > Ji/2 is 
supported by the first principle calculation Meanwhile, 
the spin susceptibility contributed from the bare band reaches 
its maximum at (n, 0) and (0,7r). Thus this relationship be- 
tween J q and the bare spin susceptibility enhances the spin 
fluctuation in the framework of RPA and correspondingly the 
superconductivity. 

We now turn to address the origin of the (n/2,n/2) spin 
fluctuation based on the fermiology picture. The FS with dop- 
ing 6 = 0.5 is plotted in Fig. 3(a). The FS sheets around X and 
its symmetric points are clearly seen. The wave vector con- 
necting the different FS sheets is just Q = (n/2, n/2). To com- 
pare the present one with the cuprates and to get an intuitive 
understanding of the interplay among the fermiology, the spin 
fluctuation, and superconductivity, we replot the FS sheets in 
the upper right BZ in Fig. 3(b). The dashed lines are the (n, 0) 
shift of the FS sheets from the upper left BZ. We can define 
the crossing points of the FS with the lines (k x = k y + ji/2) as 
the hot spots. Then the wave vector connecting the hot spots 
on different FS sheets, shown in the figure, are the wave vec- 
tor Q. The distance between the different FS sheets is close to 
Q. For the parent compound, the area of vertical shadowed 
region should equal to that of horizontal shadowed region. 
The wave vector connecting the FS sheets should be close to 
(+7r/2, +n/2). Based on the fermiology and Eq.(3), it is rather 
clear that the spin susceptibility reaches its largest value at Q 
because Sk+Q - «k is vanishingly small. For the doped sys- 
tems, the balance between the two kinds of shadow is broken, 
and the hot spots shift to the BZ boundary or diagonal direc- 
tion. The features of the FS topology and doping evolution 
seen from Fig. 3(b) are quite similar to the FS of cuprates in 
the whole BZ, except that the wave vectors connecting the hot 
spots are (+n/2, +n/2) for this compound by noting that the 
corresponding vectors are (±n, ±n) for cuprates. In Fig. 3(b), 
only two FS sheets actually exist. For cuprates four FS sheets 
exist. As a result, the electron filling for the parent compound 
is only half to that for cuprates, namely, only 0.5 for the parent 
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FIG. 4: (Color online) (a) The gap function from the self-consistent 
calculation, (b) The gap along the FS sheets with the dashed and 
solid lines are around the X and Y points, respectively, with 6 denoted 
in panel (a). 



KFe 2 Se 2 . 

Besides the fermiology, it is insightful to compare the 
present material with cuprates through the Heisenberg spin 
interaction term Hj. For cuprates, only the nearest-neighbor 
interaction needs to be considered with J\ > 0, leading to the 
AF ground state. When the AF ground state is suppressed, the 
spin excitation occurs at the wave vector near AF wave vec- 
tor (7T, 7r). For the present material, the parent compound is 
modeled by the /j — J% — J3 Hamiltonian J20ll . Here J] < 
and J2, J3 can be comparable to J\, which is quite different 
from that of cupates. The ground state is frustrated and the 
combined effect from all 7,-term leads to the block AF state, 
with one block consisting of 2 x 2 sites. The period along x 
ory direction is twice of that of cuprates. Thus when the spin 
order is suppressed, the wave vector of the spin excitation is 
half of that in cuprates, namely, (n/2,n/2). 

At this stage, we look into the pairing symmetry of the 
SC gap, which can be evaluated from Eq.(5), namely, the 
temperature T is decreased until the maximum eigenvalue 
A = 1 is obtained, with the SC gap being the eigenvector 
for the maximum eigenvalue. The intensity plot of the SC 
gap is displayed in Fig. 4(a). Obviously, the pairing sym- 
metry is of af A 2_,,:-wave. We have checked numerically that 
this result is rather robust to the reasonable parameter change. 
The SC gap along the FS for the of-wave pairing symmetry 
A = Ao/2(cosfc v - co&ky) is plotted in Fig.4(b). As is seen, 
the SC gap along one FS sheet is nearly isotropic with the 
difference less than 15%, being consistent with the ARPES 



experiments lfl4T - [T7ll . In addition, it was reported in ARPES 
experiments that the gap is small or vanish for the smaller 
pockets around the T point HUGH], being also qualitatively 
consistent with the d x i .,2-wave. However, it was indicated in 
Refs. lfl5i [I7I1 that there exists a larger FS pockets around F 
with the gap being nearly isotropic. This result seems to con- 
tradict with the (i-wave symmetry. While we think that the 
larger FS around F is actually due to the band folding effect: 
dominant features of the pocket and the gap may be folded 
from those around M point. If this is the case, our results for 
the gap symmetry agree qualitatively with the ARPES exper- 
iments. 

The origin of the d-wave superconductivity can be under- 
stood based on the spin fluctuation picture and FS topology. 
The pairing potential V contributed by the spin fluctuation 
is largest at the wave vector Q = (±n/2, ±n/2). The factor 
tanh(j6ei t -/2)/2ek' in Eq.(5) is positive for any and largest 
at £k' = 0, which means that the pairing near the FS is im- 
portant. As a result, the gap function at or near the FS should 
satisfy the condition Ak = -Ak+Q. As shown in Fig. 3, if k 
belongs to one sheet of FS, then k + Q should be near the 
other neighboring sheet of FS. For the d x i_ v 2 -symmetry, as dis- 
played in Fig. 4(b), the SC gaps have the same magnitudes and 
different signs along the two neighboring FS sheets. The con- 
dition Ak = -Ak+Q is satisfied approximately. In this sense, 
we give an intuitive understanding of the pairing symmetry in 
this material. 

In summary, we have proposed a single-band model to de- 
scribe the K A Fe2Se2 material. The band structure has been 
found to be qualitative consistent with the five-orbital model 
at low energies. The spin excitation with the wave vector 
(+7r/2, +7t/2) has been revealed. Our theoretical results can 
well be understood based on the dierct fermiology analysis. 
In particular, we have presented a sound and coherent picture 
for the interplay of the SDW and superconductivity in this 
compound and suggested that the essential physics is some- 
how similar to that in other high-T c SC materials. Finally, the 
rather robust d^—f. pairing symmetry of superconductivity has 
been revealed. 
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